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ABSTRACT 

We have conducted an optical and infrared imaging in the neighbourhoods of 4 triplets 
of quasars. R, z', J and K s images were obtained with MOSAIC II and ISPI at Cerro 
Tololo Interamerican Observatory. Accurate relative photometry and astrometry were 
obtained from these images for subsequent use in deriving photometric redshifts. We 
analyzed the homogeneity and depth of the photometric catalog by comparing with 
results coming from the literature. The good agreement shows that our magnitudes 
are reliable to study large scale structure reaching limiting magnitudes of R = 24.5, z' 
= 22.5, J = 20.5 and K s — 19.0. With this catalog we can study the neighbourhoods 
of the triplets of quasars searching for galaxy overdensities such as groups and galaxy 
clusters. 
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1 INTRODUCTION 

Our present understanding of galaxy formation and evo- 
lution can be improved significantly through studies of 
quasars. In the unified scenario, merger-driven processes in- 
clude the origin of a quasar by self-regula ting, supermassive, 
black-hole growth (|Hopkins et alj l2005h . Thus, most mas- 
sive spheroids should host an optically-bright quasar during 
a brief period of its evolutionary history. 

There is well-established evidence indicating that sev- 
eral galaxy properties depend on the environment where 
the galaxies formed and evolved. Several processes - such 
as stimulated or truncated star formation, tidal stripping, 
merging, etc. - can determine the nature of galaxy popu- 
lations and dynamics. Also to be considered is the possi- 
bility that AGN feedback processes could induce significant 
changes in t he evolution of gala xies near to the quasar (see 
for instance ICroton et alj (|2006l )b By studying quasar envi- 
ronments we may deepen our understanding of the relation 
of quasar fueling and galaxy interactions and merging; and 
as a consequence, the formation of bright galaxies. 

The availability of large quasar samples has increased 
significantly in recent years, mainly due to the advent 
of l arge surveys such as the 2dF Galaxy Redshift Sur- 
vev (IColless et all 200lh and the Sloan Digital Sky Survey 
l|Stoughton et alj|2002h . These samples provide an opportu- 
nity for more robust statistical studies of the quasar phe- 
nomenon and its link with galaxy formation. 



At low redshifts (z < 0.3), existing cross-correlation 
analysis of quasar environments show that quasars in- 
habit enviro n ments similar to thos e of n ormal galaxies 
jSmith et all i|l995l ). IColdwell etafl (|2002r i. However, on 
small scales (projected distance, r p < lh -1 Mpc; and radial 
velocity difference, AV = 500km s -1 ), the quasar environ- 
ment is overpopulated by blue disk galaxies having a strong 
star-formati on rate with respect to t ypica l galaxy neigh- 
bourhoods (|Coldwell fc Lambasll2003l, I2OO6I ) . Thi s resul t is 
in agreement with those of ISochting et alj l|2002l . |2004| ) in 
the sense that low-redshift quasars follow the Large Scale 
Structure traced by galaxy clusters but they are not located 
near the center of clusters. These quasars are mainly found 
in the periphery of such structures or between two, possibly 
merging, clusters. 

At higher reds hifts, quasars are oft e n associated with 
rich environments (|Hall fc Greenl Il998l ; iDiorgovskil Il999ft 
where the comoving density of galaxies is higher than that 
expected for the general field. This should be interpreted 
as the core of future rich clusters. Some results at z ~ 1 
suggest that quasar s reside on the peripherie s of clusters or 
in cluster m ergers llHaines et al. I l200ll . |2004 iTanaka et al.l 
2000. l200ll) . lTanaka et al l <|200lT >. for example, investigated 
a group of 5 quasars tracing a ~ 10 Mpc structure of 4-5 
clusters, with only a single radio-loud quasar appearing to 
be directly associated with one of the clusters. 

In general, the existing studies of quasar environments 
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suggest their association with forming structures, marking 
merging clusters and filaments. Consequently, groupings of 
quasars can be expected to trace regions of extraordinary 
activity. Pairs of quasars have been strongly studied and 
they are used to identify cluster of galaxies at high redshifts. 
Zhd anov fc Surdeil l|200ll ) suggested that quasar pairs can be 
used as tracers of high redshi ft large-scale structures. More- 
over, iDiorgovski et al.l (|2003h found a quasar pair at z ~ 5 
associated with large scale struct ure, possibly a pro toclus- 
ter at that redshift. More recently. iBoris et al.l l|2007l ) found 
that quasar pairs at z ~ 1 are excellent tracers of high den- 
sity environments suggesting that they can be used to find 
galaxy clusters. 

Triplets of quasars in a relative small volume are ex- 
tremely rare and until now, no cluster-scale triplets have 
been published or investigated. A very compact quasar as- 
sociation (~ 50 Kpc s eparation) was recently reported by 
IDiorgovski et all (|2007l ) suggesting that this could be a com- 
pact galaxy group in a merging process. We propose to use 
triplets of quasars as "lighthouses" to find and study clus- 
ters of galaxies that favour such rare events. With this aim 
in mind, we performed an optical and infrared photometric 
study of the environments of triplets of quasars at differ- 
ent redshifts on the scale of galaxy clusters. The images 
allow us to observe galaxies in the early Universe and they 
are used to investigate the environments of quasar triplets 
We also study the group and cluster morphology, and the 
photometric properties of individual galaxies in the neigh- 
bourhood of each triplet. We have already found that low-z 
triplets of Seyfert 1 gal axies are associated w ith extremely 
rich clusters of galaxies i|Sochting et alj|2007jl If the quasar 
triggering mechanism does not evolve strongly over redshift 
we can expect to find some of the richest galaxy clusters ever 
discovered at z ~ 1. A negative result will be clear evidence 
of the evolution of the quasar formation mechanism. 

This paper is organized as follows: In section 2 we de- 
scribe the sample selection and section 3 shows the optical 
and infrared observations and data reduction. Section 4 de- 
scribes the photometric galaxy catalog, and in section 5 we 
provide a summary of the main results. We adopt the latest 
cosmology with Qm = 0.3, Q\ = 0.7 and a Hubble constant, 
H = 70 km s" 1 Mpc" 1 . 



2 THE SAMPLE 

The quasar sample is defined by quasars which have at 
least one emission line with a full width at half maxi- 
mum larger than 1000 km s -1 , luminosities brighter than 
Mi - -23, PSF i magnitude < 19.1 and highly reliable red- 
shifts (Schneider et al. 2003). To analyze the regions in the 
neighbourhoods of groups of quasars, we identified different 
active systems using t he friend-of-friend, three- dimensional, 
percolation algorithm (|Huchra fe Gellen l|l982h ). This tech- 
nique was applied to the sample of 54000 quasars taken from 
SDSS-DR4 which satisfy the criteria defined above. Taking 
into account galaxies in clusters and their surroundings, and 
the low spatial density of quasars (about 15 quasars per 
square degree), we finally defined as candidates for triplets 
of quasars those systems with percolation longitudes smaller 
than 2 hT 1 Mpc and 2000 km s _1 . These limits are about 
twice the typical region for rich galaxy clusters. We found 



293 pairs of quasars, only 7 triplets and no systems with 
more members. To increase t he sample of triplets, we ap - 
plied the same criteria to the IVeron- Cettv fc Veronl (|2003D 
catalogue. 

Figure [1] shows the redshift distribution of the fraction 
of the quasar number normalized to the total quasar sam- 
ple, and that corresponding to pairs and triplets (systems) 
of quasars found with the adopted selection criteria. The 
redshifts of each of the twelve triplets in the total sample 
of triplets are marked at the top of the figure. For systems 
of quasars, there are two peaks at z ~ 0.2 and z ~ 1.8, 
which are more concentrated than the total quasar sample. 
In spite of the uncertainties due to the small amount of data 
for systems of quasars, the presence of these two peaks may 
be interpreted in terms of differences in t he growth of self- 
regul ated black holes for different masses (|Di Matteo et all 
120051 ). They proposed a fast evolutionary growth for the 
massive black holes formed in the early Universe, which ex- 
hausted the surrounding gas due to a higher accretion rate. 
On the other hand, the lifetime of the active phase increases 
for black holes of smaller masses implying an excess of low 
luminosity quasars in the local Universe. 

The total sample consists of 12 triplets of quasars: 3 
triplets with z < 0.2 and 9 within the range of 0.9 < z 
< 2.5. As mentioned above, pairs of quasars - as a class - 
have been extensively studied and so we will concentrate our 
efforts here only on triplets of quasars, which are very rare. 

Our first subsample of triplets o f quasars, those nearb y 
triplets with z < 0.2 is discussed in ISochting et afl l|2007h . 
Using the photometric and spectroscopic catalogs of the 
SDSS-DR4, we found a strong association of z < 0.2 triplets 
with the richest central parts of superclusters. In all cases 
the individual members of the triplets have been found on 
the periphery of the main galaxy concentrations. 

For some of our second subsample, those triplets at 
higher redshifts, we performed multicolor photometry of the 
regions near to these triplets in order to locate associated 
galaxy groups and clusters. 



3 OBSERVATIONS AND DATA REDUCTION 

In this paper, we concentrate on four high-redshift triplets 
from the sample. We obtained multicolor photometry of 
their neighbourhoods with the CTIO 4m telescope using 
two instruments: MOSAIC II for the optical and ISPI for 
the near IR (Proposal Ids. 06A-0146 and 06B-0117). Ta- 
ble 1 shows the observed triplets, including the triplet and 
quasar identifications, coordinates, redshifts and absolute V 
magnitudes, respectively. 

3.1 Optical imaging 

Triplets 5 and 6 were observed in the R and z' bands using 
MOSAIC II during 5 half nights from February 2 nd through 
6 , 2006. This optical instrument is an array of 8 2048 x 
4096 SITe CCDs with a gain of 2 e-/ADU and a scale of 0.27 
arcsec/pixel, giving a total field of view of 36 x 36 arcmin 2 . 
In the R band, we made eight exposures of 600s for Triplet 
5 and eleven exposures of 600s for Triplet 6. Moreover, four 
exposures of 600s were taken in z' band for both triplets. 
Table 2 shows the triplet Ids together with a summary of 
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the optical observations, in R and z', in columns 2 and 3, 
respectively. 

3.1.1 Photometry 

The images we r e reduced using the SMSN p i peline 
jRest et al.1 120051 ; iGarg et all 120071 ; iMiknaitis et~atl I2007T I 
originally created for the SuperMACHO and ESSENCE 
projects. The SMSN pipeline reduces the raw images using 
the IRAF Q mosaic routines XTALK and MSCCMATCH, 
as well as its own routines in C, perl and python. As part of 
this reduction, the images are crosstalk corrected, astromet- 
rically calibrated, and split into the 16 separate amplifiers. 
The images are then bias subtracted and flat fielded. The 
z-band images were also defringed. 

We flux calibrated the R a nd z' image s using our 
observations of several lLandoltl (|l992f) and ISmith et al.l 
(2002) standard stars, respectively. In our images, 
point-source magnitudes were obtained using SExtractor 
IjBertin fc ArnoutsHl996T ) with an aperture diameter of 6", 
which proved to be adequate by monitoring the growth curve 
of all the measured stars. For the final calibration, we in- 
clude the zero point and extinction term. The colour term 
was negligible compared to m a gnitu de errors and it was not 
included (as in lGawiser et al.l ((2006)). 

For R magnitudes, the zero point is in agreement with 
the published value for the instrument. In order to better 
analyze the precision of our calibration, we also compared 
the magnitudes of the stars in our fields obtained with SEx- 
tractor (MAG_BEST) with the PSF magnitudes computed by 
SDSS. Before performing the comparison we transformed 
SDSS r' magnitudes to our system using Fadda et al. (2004) 
transformations. Figure[2]shows in the upper panel this com- 
parison of R magnitudes showing the good agreement with 
those of the SDSS. The mean difference is A R = 0.162 ± 
0.139. 

In the z'-band the comparison with SDSS was made 
without any transformation because the filters are identical. 
The bottom panel of Figure [2] shows this comparison. The 
mean difference is A z' = 0.14 ± 0.08 mag. This offset is 
consistent with the estimated uncertainty of the overall zero- 
point, mainly due to the very small number of standard stars 
used in this band. Our final z' magnitudes were corrected by 
this offset. 

We applied the offsets in both bands to correct our final 
magnitudes. Note also that all magnitudes are in the Vega 
system. 

3.1.2 Astrometric precision 

In order to assess the accuracy of the astrometric calibration, 
we compared the position of the stars in the observed fields 
with those of the USNO-A2.0 Catalog (Monet et al. 1998). 
We used only unsaturated stars from our images. Figure [3] 
shows only the results of the comparison for the z'-band. 
The mean offset in the positions between both catalogs is 
negligible: Aa = -0.057 ±0.281 arcsec and AS = -0.097 ± 

1 Image Reduction and Analysis Facility (IRAF), a software sys- 
tem distributed by the National Optical Astronomy Observatories 
(NOAO) 



0.302 arcsec for R-band and Aa = —0.110 ± 0.236 arcsec 
and AS = -0.074 ± 0.302 arcsec for z'-band. The scatter 
in the coordinates is found to be typically smaller than 0.3 
arcsec and it is consistent with the astrometric accuracy of 
the USNO-A2.0 catalog. We also compared the positions of 
the stars with SDSS sources and we found systematic offsets 
of 0.218 ± 0.096 and 0.434 ± 0.107 in right ascension and 
declination, respectively - similar to those found by Fadda 
et al. (2004). 

3.2 Infrared imaging 

The triplets of quasars were also observed in the infrared us- 
ing ISPI in two runs during 2006: 5 half nights from April 2 nd 
to 6 th and 3 complete nights - October 7 th to 9 th . ISPI (In- 
frared Side Port Imager) is a 1-2.4 [im imager with an 2048 
x 2048 HgCdTe array with an scale of 0. 3 arcsec /pix, giving 
a fiel d of view of 10.5 x 10.5 arcmin 2 (|van der Bliek et al.l 
[2004 ) . The gain is 4.25 e-/ADU and the dark current is very 
low - about 0.1 e-/sec for long exposures. However, this ar- 
ray has a number of bad pixels which are located along a 
few rows and columns in the upper right corner. 

All four triplets were observed in the K s band but only 
Triplets 6 and 7 were observed in the J band. We used a cir- 
cular dithering pattern of 15 different positions to observe 
the triplet fields. For each image in the sequence, the expo- 
sure time was 130s x 1 coadds and 20s x 5 coadds for J 
and K s , respectively. Table 2 includes a summary of the IR 
observations in columns 2 to 5. 

3.2.1 Photometry 

For the reduction, we used standard tools in IRAF and some 
tasks from the CTIO Infrared reduction package (cirred) 
written by R. Blum for specific ISPI processing. The basic 
data reduction consisted in creating the master flatfield im- 
age from a median of lamps on - lamps off dome flats and 
the construction of the bad pixel mask. The final master flat 
was corrected by illumination effects and normalized by the 
average intensity. Median sky images were created for each 
dithered sequence. The scientific images were fiatfielded, bad 
pixel corrected, and sky subtracted. T he astrometr y was per- 
formed with the package WCSTools0 (|Min]j2006h ■ Then for 
each sequence the images were weighted combined with indi- 
vidual weight map images using SWarp written by E. Bertin 
B 

To calibrate the magnitudes, we observed some stan- 
dard stars from Persson et al. (1998). We also used those 
stars in our fields belon ging to the 2MASS point source cat- 
alog (|Cutri et al.ll2003T ) to improve precision. As for the op- 
tical magnitudes, the stellar instrumental magnitudes were 
obtained with an aperture diameter setting of 6". For the 
final calibration, again we did not include the color term. 

To verify the photometric precision we compared our 
magnitudes of detected objects (those with 5 a above the 
sky background) with the corresponding, uncontaminated, 
point sources from the 2MASS catalogue. Figure [4] presents 
comparisons of the magnitude differences for the stars in 

2 Available at ftp://cfa-ftp.harvard.edu/pub/gsc/WCSTools/ 

3 http://terapix.iap.fr/soft/swarp 
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common with 2MASS in J and K s bands, respectively; the 
agreement is excellent. In the J band, for objects brighter 
than 16 mag, A J = -0.007 ± 0.153 mag. In the K s band, 
the comparisons look similar; AK S — — 0.066±0.109 mag for 
objects brighter than 14.5 mag. Finally, in both bands, we 
can see from Figure [4] a progressive increase in the scatter 
at fainter magnitudes due to larger errors in the 2MASS 
catalogue at these magnitudes. 

3.2.2 Astrometric precision 

In order to verify the precision of our astrometry, we also 
checked our measured positions for 2MASS objects in our 
fields against the 2MASS catalogue. Figure [7] shows the dif- 
ferences in both coordinates, right ascension and declination, 
between the two sets of data in K s band. The average differ- 
ences are Aa = 0.096 ±0.396 arcsec and AS = 0.006 ±0.347 
arcsec for J band for objects brighter than 16 mag; and 
Aa = 0.038 ± 0.313 arcsec and AS = 0.018 ± 0.371 arcsec 
for K s band for objects brighter than 14.5 mag, an overall 
astrometric precision sufficient for our purposes. 



4 SOURCE DETECTION AND CATALOGUE 

Our photometric catalog in the four bands was constructed 
using SExtractor. An object was considered detected if the 
flux had an excess of 1.5 times the local background noise 
level over at least 5 connected pixels. 

The catalog includes magnitudes in apertures of 2, 3 
and 4 arcsec, total magnitudes, star-galaxy separation, ellip- 
ticity, etc. The star-galaxy classification represented by the 
stellarity index, CLASS_STAR allows us to decontaminate the 
galaxy catalog. Figure [6] shows this parameter as a function 
of R magnitudes. The results are similar for z' band. From 
the Figure we can see clearly two object sequences: stars de- 
fined with CLASS_STAR ~ 1 and galaxies with CLASS.STAR ~ 
0. At fainter magnitudes, Bertin & Arnouts (1996) showed 
that the stellarity index is strongly seeing dependent. Tak- 
ing this into account, we defined as galaxies those objects 
with CLASS_STAR < 0.9 and CLASS.STAR < 0.75 for R and z' 
bands, respectively. 

The IR photometric catalog was constructed using SEx- 
tractor in double-image mode with the K s image as a refer- 
ence on account of its better quality. 

Figure [7] shows the IR colour-magnitude diagram and 
how the J — K a colour alone provides a reasonable means 
of performing the star-galaxy separation. The horizontal se- 
quence at J — K s < 1 corresponds to the stellar sequence 
of g alactic stars with ty pes later than G5 and earlier than 
K5 (|Finlator et al J 12000). In the Figure we also plot point 
sources detected with CLASS_STAR > 0.95 (star symbols). As 
can be seen most of these objects are located within the stel- 
lar locus box. But those point sources with redder colours 
(J — K s >1 mag and K s > 18) are too faint and the star- 
galaxy separation becomes unreliable. In spite of this, we 
prefer to include them in the final catalog (also using star 
symbols) . 

To reinforce our star-galaxy separation we also use the 
half-light radius parameter, ri/ 2 . This parameter measures 
the radius that encloses 50% of the object's total flux. It 
is clearly independent of magnitude and depends only of 



the image seeing. In Figure [8] we show the half-light radius 
vs K s magnitude. Again, as in Figure [7] point sources with 
CLASS_STAR > 0.95 are plotted with star symbols. As men- 
tioned in Sub-Section Infrared imaging, there are some bad 
pixels located along a few rows and columns and yet they 
were not rejected during the reduction process. These ob- 
jects have r 1 / 2 < 1 in the Figure and they were discarded 
in the final catalog. Also, the stellar locus is clearly visible 
at rx/2 ~ 2 and K a < 18. Therefore, we define as galaxies 
those objects with J — K s > 1 and K s < 18 mag. 

We created the final photometric catalog discarding 
spurious objects and stars. 

4.1 Results. Number counts 

The final catalog in all bands was obtained with total mag- 
nitudes, and no need to apply an aperture correction. Our 
fields are located at | b |> 40°, so the Galactic extinction 
effect is generally small. Only the R magnitudes were cor- 
rected for galactic extinction using the E(B — V) values pro- 
vided in the NASA / IPAC Extragalactic Database (NED) 
Q, which are based on the E(B — V ) values from the extinc- 
tion maps of ISchlegel et all (I1998T ). 

The number counts as a function of magnitude allow us 
to analyze the homogeneity and depth of the catalog and to 
compare it with results obtained from the literature. 

Table 3 shows the main result, the number of galaxies 
per square degree per magnitude in the four bands: R, z', J 
and K s and their respective errors. The Figures EJl 1 10 1 [TT1 and 
1121 show the galaxy number counts, in the four bands. We 
compared our results with galaxy number counts from the 
literature, mainly: MacDonald et al. (2004) and Metcalfe et 
al. (2001) in the R band; Kashikawa et al. (2004) and Capak 
et al. (2004) in the z' band; Olsen et al. (2006), Saracco 
et al. (1999) and Iovino et al. (2005) for the J band; and 
finally Best et al. (2003), Olsen et al. (2006), Best (2000) and 
Bornancini et al. (2004) for K s band. From these Figures it 
is clear that we reach magnitudes of R = 24.5 mag, z' = 22.5 
mag, J = 20.5 mag and K s = 19.0 mag. Within these limiting 
magnitudes, we find good agreement with the galaxy counts 
coming from other authors in all bands. 



5 SUMMARY 

We selected a sample of triplets of quasars as those triple 
systems with percolation longitudes smaller than 2 ft -1 Mpc 
and 2000 km s _1 . The goal is to use them as "lighthouses" 
to find richest clusters. For this purpose we conducted a 
photometric study of the neighbourhoods of 4 triplets of 
quasars in the redshift range of 0.9 to 1.7. Images in R, 
z', J and K s bands were obtained with MOSAIC II and 
ISPI at CTIO. Magnitudes and coordinates were found to 
be in excellent agreement with the USNO-A2.0 and 2MASS 
catalogs. We analyzed the homogeneity and depth of the 
photometric catalog by comparing the number counts with 
results coming from the literature. The excellent agreement 
shows that our magnitudes are reliable, down to limiting 

4 http://nedwww.ipac.caltech.edu/ - the NASA-IPAC Extra- 
galactic Database 
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magnitudes of R = 24.5, z' = 22.5, J = 20.5 and K s = 19.0. 
These results allow us to study large- scale structure and in 
a forthcoming paper we will deal with the statistical analysis 
of this data. 
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Table 1. The Observed Sample of Triplets of Quasars 



Triplet 
Quasar Id. 


RA(J2000) 


DEC(J2000) 


z 








TRIPLET 4 




















2QZ J015646-2828 


1 


56 


46 


6 


-28:28:0 





919 


-24 


5 


2QZ J015647-2831 


1 


56 


47 


9 


-28:31:43 





919 


-24 


4 


Zi Zj J IJ J_ tJUtJU iiOZiO 


1 


56 


56 





- .^o .oo 


o 


925 


-24 


3 


TRTPLFT 5 




















SDSS J10023+0155 


10 


02 


19 


5 


01:55:37 


1 


509 


-24 


8 


SDSS J10025+0150 


10 


02 


34 


3 


01:50:10 


1 


506 


-25 


7 


SDSS J10026+0159 


10 


02 


36 


7 


01:59:48 


1 


516 


-24 


7 


TRIPLET 6 




















2QZ J120914+0035 


12 


09 


14 


9 


00:35:51 


1 


316 


-25 


8 


2QZ J120919+0029 


12 


09 


19 


6 


00:29:26 


1 


319 


-24 


3 


2QZ J120922+0026 


12 


09 


22 


4 


00:26:46 


1 


322 


-24 


2 


TRIPLET 7 




















2QZ J235405-2845 


23 


54 


05 


3 


-28:45:06 


1 


673 


-25 


2 


2QZ J235430-2848 


23 


54 


30 


2 


-28:48:41 


1 


670 


-24 


6 


2QZ J235414-2848 


23 


54 


14 


2 


-28:48:10 


1 


679 


-25 


8 



Table 2. Observing Logs 



Triplet Id. R z' J K, 



TRIPLET 4 

TRIPLET 5 8 x 600s 4 x 600s 
TRIPLET 6 11 x 600s 4 x 600s 
TRIPLET 7 



7 x 15 x 5 x 20s 

8 x 15 x 5 x 20s 
15 x 15 x 130s 8 x 15 x 5 x 20s 
12 x 15 x 130s 7 x 15 x 5 x 20s 



Table 3. Number Counts in N Deg -1 0.5 Mag" 1 



Magnitude 


R 


z' 


J 


K s 


14.75 








126 ±63 


15.25 








347 ±104 


15.75 








663 ±144 


16.25 






410 ±113 


1359 ±207 


16.75 






726 ±151 


2117 ±258 


17.25 




151 ±21 


979 ±175 


4013 ±356 


17.75 




281 ±29 


2180 ±262 


5466 ±415 


18.25 




446 ±36 


2054 ±254 


10110 ±565 


18.75 




814 ±49 


4013 ±356 


13680 ±657 


19.25 


357 ±32 


1489 ±66 


6098 ±438 


15580 ±701 


19.75 


804 ±48 


2269 ±81 


7520 ±487 


15580 ±701 


20.25 


1416 ±64 


3827 ±105 


11630 ±606 


13020 ±641 


20.75 


2112 ±78 


5946 ±131 


12950 ±639 


5814 ±428 


21.25 


2945 ±92 


8487 ±157 


13110 ±643 




21.25 


2945 ±92 


8487 ±157 


9826 ±557 




21.75 


4599 ±116 


14490 ±205 


6667 ±459 




22.25 


6754 ±140 


18820 ±233 






22.75 


10340 ±173 


17250 ±224 






23.25 


16160 ±217 


10180 ±172 






23.75 


26070 ±275 








24.25 


37260 ±329 








24.75 


40550 ±343 








25.25 


11050 ±179 
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Figure 1. Relative redshift distribution of the total quasar sam- 
ple (dashed line) , and systems of quasars represented by pairs and 
triplets (solid line). The redshifts of the 12 triplets in our sample 
are marked at the top of the figure. 



Figure 3. Astrometric differences in coordinates between the 
matched MOSAIC II and USNOA2.0 stars in the z'-band. 





Figure 2. Magnitude comparisons for optical bands between 
stars in our observed fields and those stars in common with SDSS. 
Upper panel shows the R-band magnitude differences before ap- 
plying Fadda ct al (2004) transformation to SDSS r' magnitudes. 
Bottom panel shows the z'-band magnitude differences. 




Figure 4. Magnitude comparisons for IR bands between stars in 
common with 2MASS. Upper and bottom panels show magnitude 
differences vs our magnitudes for J and Ks bands, respectively. 
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Figure 5. Differences in right ascension and declination of stars 
in common with the 2MASS catalog in the K s band. 




Figure 6. The star-galaxy separation versus total magnitudes for 
R band. 
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Figure 7. IR Color-magnitude diagram to discriminate among 
stars and galaxies. Star symbols correspond to those objects with 
CLASS.STAR > 0.95. 
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Figure 8. Half-light radius vs K s magnitude. As in previous 
Figure, star symbols correspond to those objects with CLASS.STAR 
> 0.95. The vertical line marks our dctcctability limit. 
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Figure 9. Number counts N(R) in the R band. We also show 
number counts coming from other authors: MacDonald et al. 
(2004) and Metcalfe et al. (2001). 



Figure 11. Number counts N(J) in the J band. We also show 
number counts coming from other authors: Olsen et al. (2006), 
Saracco ct al. (1999) and Iovino ct al. (2005). 
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Figure 10. Number counts N(z') in the z' band. We also show 
number counts coming from other authors: Kashikawa ct al. 
(2004) and Capak et al. (2004). 



Figure 12. Number counts N(K S ) in the K s band. We also show 
number counts coming from other authors: Best et al. (2003), 
Olsen et al. (2006), Best (2000) and Bornancini ct al. (2004). 



